Quantum-Confined Stark Effect in polar and nonpolar Wurtzite InN/GaN 
Heterostructures: Influence on Electronic Structure and Compensation by Coulomb 

Attraction 



Stefan BartheL/'F] Kolja Schuh, 1 Oliver Marquardt, 2 ' 3 Tilmann 
Hickel, 3 Jorg Neugebauer, 3 Frank Jahnke, 1 and Gerd Czycholl 1 

1 Institute for Theoretical Physics, University of Bremen, Germany 
2 now at: Paul-Drude-Institut fur Festkorperelektronik, Berlin, Germany 
3 Max-Planck Institut fur Eisenforschung, Diisseldorf, Germany 
(Dated: January 14, 2013) 

In this paper we systematically analyze the electronic structures of polar and nonpolar wurtzite- 
InN/GaN quantum dots and their modification due to the quantum-confined Stark effect caused by 
intrinsic fields. This is achieved by combining continuum elasticity theory with an empirical tight 
binding model to describe the elastic and single-particle electronic properties in these nitride systems. 
Based on these results, a many-body treatment is used to determine optical absorption spectra. The 
efficiency of optical transitions depends on the interplay between the Coulomb interaction and the 
quantum-confined Stark effect. We introduce an effective confinement potential which represents the 
electronic structure under the influence of the intrinsic polarization fields and calculate the needed 
strength of Coulomb interaction to diminish the separation of electrons and holes. 



I. INTRODUCTION 

The quantum-confined Stark effect (QCSE) is the main 
reason for a reduced optical recombination efficiency in 
wurtzite nitride heterostructures grown along a polar 
axis. For instance, quantum dots will exhibit a charge 
carrier separation, see e.g. Ref . (fTj) . since strong intrinsic 
electrostatic fields caused by piezo- and pyroelectric con- 
tributions drastically modify the confinement potential 
for electrons and holes. Motivated by the fact that the 
alloyed nitride systems can principally cover the whole 
visible range of emission wavelengths, one possible idea 
to overcome the drawback of the charge separation is to 
perform growth along a nonpolar facet, Refs.(2-4). As a 
result, the spontaneous polarization part of the intrinsic 
fields can be reduced, which should lead to weaker spatial 
separation of electron and hole wave functions and thus 
to more efficient light emission. Driven by the large po- 
tential as an application in optoelectronic devices, see 
Ref. ([5]), we theoretically analyze the electronic struc- 
ture of polar (c-plane) and nonpolar (m-plane) InN/GaN 
quantum dots (QDs) under the influence of the persistent 
intrinsic electrostatic fields. To establish a more general 
understanding, we not only discuss ground state prop- 
erties which have been studied together with the intrin- 
sic fields in Refs.([TJ and[7|, but also excited states 
of the system. These excited states play a major role if 
one includes Coulomb interaction and goes from a single- 
particle description to a many-body theory, which en- 
ables to study the competition between intrinsic fields 
and electron-hole attraction. The authors have already 
shown in Ref. (|36|) . that the Coulomb interaction does en- 
hance the optical dipole strength in pure lenshaped non- 
polar InN/GaN QDs. Recently, Schulz et ai. have stud- 
ied nonpolar GaN/AIN QDs by taking into account the 
precise QD geometry and using a self-consistent Hartree 
approximation to account for Coulomb effects in Ref. 



Their calculations for the increased electron and hole 
wave function spatial overlap is in agreement with PL 
measurements in Ref. ((3]) and supports the importance of 
geometry and Coulomb interaction for accurate model- 
ing of the optical properties of a-plane GaN QDs. This 
manuscript is intended to provide theoretical and compu- 
tational details of the approach the authors have used in 
Ref. (f36|) . in which parts of the results have already been 
published. Moreover the present manuscript discusses 
additional findings. 

In the first part of the paper we present the geometry, 
shape and composition of the systems. Subsequently, we 
provide the theoretical background on the methodology 
applied in our approach and then end with the discussion 
of the results. To begin with, the approach to approxi- 
mately determine the number of bound states within QDs 
is presented and it is explained how a calculation of an 
effective confinement potential is seperately possible for 
electrons and holes. The two-dimensional, single-particle 
effective potential takes into account the geometry, in- 
trinsic fields and confinement due to the material sys- 
tem. Based on this quantity, an understanding of the 
wave functions's localization will be gained. In Sec. |IV| 
we consider Coulomb interaction between the carriers in 
order to study the modification of carrier localization due 
to electron-hole attraction. Furthermore, we determine 
the critical interaction strength where the Coulomb in- 
teraction becomes dominant over the intrinsic elctrostatic 
fields. 



II. GEOMETRY 

In this paper we present a theoretical analysis of pure 
lens-shaped InN/GaN QDs on a wetting-layer with a 
thickness of two lattice constants in the polar (c-planc) 
and nonpolar (m-plane) growth direction. For the po- 
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lar QD, we chose a geometry according to recent experi- 
ments, Refs.([ MTTj) . In order to separate the effect of the 
built-in fields from modifications of the QD geometry, we 
use the same shape and QD dimensions for the non-polar 
QD. Note that these properties of the QDs may change if 
the heterostructure is grown in a different crystal orien- 
tation, but we did not include these geometrical effects in 
the study in order to focus on the importance of Coulomb 
interaction solely as done in Ref.(36). Nevertheless, due 
to the different crystal orientations and the discrete spa- 
tial lattice in our theory, the exact boundaries of the 
simulated nanostructure still differ within 10% between 
both geometries. The dimensions of the polar (nonpo- 
lar) QD are Ax sa 7.74 nm (Ax m 2.76 nm), Ay w 7.66 
nm (Ay sa 7.66 nm) and Az sa 3.11 nm (Az sa 7.26 
nm). In addition, the wetting-layer without the QD is 
also modelled in a seperate treatment for comparison. 
For the polar orientation we use a hexagonal supercell in 
combination with periodic boundary conditions to close 
the system and for the nonpolar orientation a cuboid 
one respectively. Throughout the paper we neglect spin- 
orbit splitting, as it is weak for InN with A so = 5 meV, 
Ref.CCi). 



already been performed in Ref.(19) for the hexagonal 
Bravais lattice by using the sp 3 -basis and taking into ac- 
count matrix elements up to second nearest neighbors in 
combination with input parameters known from GoWq 
calculations for wurtzite-nitride systems, for further de- 
tails see Ref. J20 - 23) . The dispersion e(k) throughout the 
whole Brillouin zone obtained within this model provides 
a quite reasonable description of the bulk electronic prop- 
erties, since the conduction and three valence bands are 
pinned at all critical point energies except the X-point to 
ab-initio results while in addition the anisotropic effective 
masses at the T-point are included as well as the Kane- 
parameter. More details for the quality of the analyti- 
cal parametrization can be found in Ref. (fl"9|) .We proceed 
with the respective equations in a compact form. 

In order to model a bulk system we start by defin- 
ing an eigenfunction |k) Bulk as a linear combination with 
unknown coefficients c Q (k) of Bloch sums |k, a) given by 
the Fourier sum of the ETBM basis |R, a): 



|k) Bulk = £ c «( k )^E e * kR l R 

a * iV R 



(1) 



|k,a) 



III. THEORY 

A. Electronic properties 

For the calculation of the electronic properties we use 
an empirical tight-binding model (ETBM), because it al- 
lows for the possibility of a microscopic description of 
nanostructures like QDs and can be applied to larger sys- 
tems. This is mainly due a spatial discretization on the 
scale of lattice sites in combination with the use of a min- 
imal basis set of orbitals describing the relevant chemical 
bonds. In the past, ETBMs have been successfully used 
to describe the electronic states of bulk semiconductors 
and low-dimensional heterostructures such as quantum 
wells, nanorods and QDs., see e.g. Refs. (fl3HT7|) . In 
the following, a brief summary of the approach used in 
this paper is provided. Our starting points are Loewdin- 
orthogonalized, effective local orbitals |R, a), localized 
on the Bravais lattice at position R. with a combined 
spin and orbital index a which are Wannier-like func- 
tions w(r — R). This type of ETBM is also referred 
to as an effective bond orbital model (EBOM), since it 
lacks the full atomistic resolution, but on the other hand 
there is some freedom to choose the tight binding basis 
as only hopping matrix elements E^ 1 ^, between these 
basis functions are needed. In principle, one can con- 
serve the microscopic symmetry by extending the orbital 
basis as done in Ref- lfTB]) . In general, one has to per- 
form a numerical fit or derive analytic expressions for 
-^r. R/ m terms of known physical parameters (effective 
masses, Luttinger parameters, bandgaps, critical points' 
energies, Kane parameters, etc.) either known from ex- 
periments or ab-initio calculations. The whole task has 



The corresponding matrix eigenvalue problem is of the 
form 

ff Bulk (k) = i- V e ik ( R - R ') (R', a '\H\R, a) (2) 
N Jrr?. >• ' 



R.R' 



and provides the dispersion e(k) and coefficients c a (k) 
after exact numerical diagonalization of Eq.Q. TV de- 
notes the number of unit cells and ensures proper nor- 
malization. The reader is again referred to Ref. (fl"9|) for 
the details of this procedure. 

The next step to generalize Eq.([2) for the descrip- 
tion of a two-dimensional translational invariant system, 
i.e. a wetting-layer, is done by restricting the Fourier 
sum in Eq.([2| to the translational invariant in-plane di- 
mensions solely. This approach corresponds to mapping 
the problem to a one-dimensional column whose effec- 
tive real-space hopping matrix-elements are given by a 
partial Fourier sum, or simply speaking, we factorize the 
problem into an in-plane (_L) and out-of-plane (||) con- 
tribution. This leads to the following equation in which 
iVj_ is the number of in-plane unit cells: 



(R/l +R\ l ,a'\H\R 1 _ +Ri|,a}. 



(3) 



There are additional matrix elements between differ- 
ent sites R|| , Rj| and not only between different orbitals 
a, a' compared to Eq.|2|. Diagonalization of this ma- 
trix eigenvalue problem provides the subband dispersion 
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e(kj_) and the coefficients cr, a (kj_) of the respective 
eigenstates: 

ik ± ) W1 = E c ^n ^( k ^) 7^ E eZk±R± + R n . a ) ■ 

R||,a * - 1 Ri 

v v ' 

|k x ,R||,o> 

Setting up the Hamiltonian for a fully non- 
translational invariant system, i.e. a quantum dot, is 
now straightforward and yields: 

H Qd = (R',a'\H\K,a). (5) 

The eigenfunctions | -0) are linear combinations of the 
ETBM basis with coefficients cr q : 

i^) Qd = E c ^i R ' a )- ( 6 ) 

R,a 

At the interfaces of two different materials A/B of a 
heterostructure we use the average hopping of the two 
constituents. Furthermore an appropriate relative va- 
lence band offset AEy for the on-site matrix elements 
of the two materials has to be incorporated in order to 
correctly model a heterostructure. In this work we use 
the valence band offset recommended in Ref . ()T2| which is 
0.5 eV for the InN/GaN material system. The numerical 
diagonalization of Eq.([5]) is performed using the folded 
spectrum method of Ref . (f2"%|) . 

B. Strain and elastic properties 

The elastic properties of the QDs under consideration 
have been computed using a plane-wave based implemen- 
tation of second-order continuum elasticity theory as de- 
scribed in Ref. (25). Within the wurtzite crystal struc- 
ture, the elastic energy 

F = \J v d3rC ^( e lx + 4) + C 3 3 4 

^xx^yy xx H~ £yy ) 

+ 4C 44 (eL + 4 Z ) + 2{C n - C 12 )e 2 xy (7) 

is minimized with respect to the displacements u(r). 
Here the Cij = Cy(r) are the elastic constants of the 
wurtzite structure and V is the simulation cell volume, 
Ref. (f26|) . The strain tensor %-(r) is related to the dis- 
placements via: 

with the local intrinsic strain being defined by the lattice 
constants of the materials involved, 

o fx xx ^rcf - «(r) c rcf - c(r) 

J a(r) c(r) 



where a re f and c ro f are the reference lattice constants, 
chosen to be the bulk lattice constants of the surround- 
ing GaN. 

Once the strain tensor 6y(r) is known, the built-in po- 
larization potential can be computed, arising from a 
piezoelectric and a spontacnous polarization. Within the 
wurtzite lattice, the polarization is given by 

/ 2e 15 e xz \ ( \ 

P = 2e 15 e yz + . (fO) 

\ e 31 (e xx + e y y) + e 33 e zz J \ P sp J 

Here e^- = e^- (r) denotes the piezoelectric constants and 
P,p is the spontaneous polarization in the wurtzite struc- 
ture. From the polarization vector P = P(r), the po- 
larization potential 4>p(r) is determined by solving the 
Poisson equation 

K VMr)V</> P (r)] = Mr), (11) 

where the polarization charge density is calculated as 
Qp(r) = — VP(r). The influence of electromechanical 
coupling between strain and polarization potentials^, 
was neglected since they are expected to be small. In- 
deed, only minor influences of this effect have been found 
in InGaN-based nanostructures^. The polarization po- 
tential then enters the modelling of the electronic prop- 
erties of the QD as an additional potential contribution. 
For the calculation of the elastic and piezoelectric prop- 
erties, the lattice, elastic and piezoelectric constants from 
Ref. (fT2"[) have been employed. It is here important to note 
that the choice of different piezoelectric constants (e.g. 
those from Ref. (|32|) can induce significant modifications 
to the built-in electrostatic potential. In particular, the 
choice of the parameter ei5 that has been controversely 
discussed recently ([7]) can lead to strong modifications 
mainly for the case of the nonpolar system. 

The inclusion of strain and intrinsic fields into the cal- 
culation of electronic properties can be accomplished on 
different levels of sophistication. In principle, the hop- 
ping matrix elements E^, have to be altered due to 
the strain in a heterostructure: i) If the relaxed atomic 
positions which minimize the elastic energy of the nanos- 
tructure are known, one can use the new bond-lengths 
and angles to rescale the hopping matrix elements ac- 
cording to Harrison's law^, though it might not be valid 
for our model, ii) In addition, the widely used two-center 
decomposition 13 allows for a direct analytical dependence 
of the hopping matrix elements as a function of bond- 
lengths and angles, iii) As a special case, the hamilto- 
nian of the EBOM reproduces the hamiltonian of k • p- 
theory in the limit of k = by construction. Thus, it 
may be possible to derive an analytical parametrization 
for the -Err/ which directly includes a dependence on 
the deformation potentials^ of the strain related k • p 
HamiltoniarPS. 

A second important part arises from the intrinsic elec- 
trostatic potential energy, 

V p (r) = -e-<f> P (T), (12) 
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which is naturally included as an on-site contribution 
to the tight-binding hamiltonian. This requires to first 
map and interpolate V p (r) obtained from continuum elas- 
tic theory to the tight-binding lattice R and then after- 
wards a symmetrization step. By rotating V^(R) six- 
times according to the Cq v symmetry of the polar geom- 
etry (c-plane) and calculating a mean potential V^,(R), 
any artificial symmetry spoiling can be avoided. The 
same procedure is used for the nonpolar geometry (m- 
plane) with the difference that an inversion is used as 
a symmetry operation. In this work we only consider 
the contribution arising from V^,(R), since we assume it 
to be the more siginificant contribution compared to the 
modification of the hopping matrix elements. 

C. Interacting many-body problem 

The solution of the many-body problem is found by 
employing the full configuration interaction (FCI) frame- 
work, which is a method that has been successfully ap- 
plied to QD systems with different basis setsSBESl. A 
main advantage over single-particle approximations like 
Hartree-Fock is, that all Coulomb correlations are in- 
cluded in the approach within the chosen basis. The 
reader is referred to Ref . ([3"5"]) for a more detailed discus- 
sion of this topic. The only limiting factor is the size 
of the single-particle basis, which restricts to low car- 
rier densities. In our case, we are interested in studying 
the optical properties associated with one electron-hole 
pair, so there is only the electron-hole Coulomb interac- 
tion V ehhe present and scattering processes that change 
the number of carriers are not considered in our calcu- 
lations. Furthermore, the small electron-hole exchange 
yeheh j g ne gj ec ted. The corresponding excited electron- 
hole (two-particle) eigenstate can be written as a linear 
combination of electron-hole basis states, 

NM=E C ^>>, (13) 

where e£, and hp are electron and hole creation operators 
and |0) is the vacuum state. 

In order to visualize the results of our calculations for 
the excitonic many-body eigenf unctions \ipx): it is con- 
venient to independently provide the electron and hole 
part. Since a factorization of the density-operator 

P= E 4^|0)c^4>(0|Vea'- (14) 

a,0,a',/3' 

is not possible, we take a partial trace over the electron 
or hole part in order to define a density-matrix for an 
electron (and corresponding for the hole) as 

a, a' /3 



so that the probability density is finally given by 

(R|p e |R) = E<R|aK,o<<«'|R>> (16) 

a, a' 

with I a) = 4|0) and |/3) = h\\G). The many-body 
Hamiltonian, 

h x =j2 E ^U a +J2 E ^ h V h p ( 1? ) 

a fj 

- E (V$fafea>hlh fi , + h.c.), (18) 

aPa'0' 

contains the single-particle energies E a , E» and the 
electron-hole Coulomb matrix elements V£pp? a ,, which 
are evaluated in Fourier space, 

V$fca, = ^ E ^ H ^ I"') (A l^'> ( 19 ) 

with the system volume V, elementary charge e and per- 
mittivity k. Alternatively, they can be calculated in real- 
space using the approximations in Rcfs.(33). If one re- 
stricts | a) and \j3) to a finite subspace of electron and 
hole single-particle states, which are localized in the QD 
region, a direct diagonalization of Hx becomes possible. 
In addition, the restriction to one electron hole pair al- 
lows for using a very large single-particle tight-binding 
eigenbasis entering the FCI calculation and in particular 
all bound QD states can be included this way. 

IV. RESULTS 

A. Subbandstructures 

One of the aims of this work is the investigation of ex- 
cited states. For this purpose one has to know how many 
bound QD states exist in the given nanostructure. One 
solution is to calculate an energetic cutoff which sepa- 
rates all bound single-particle QD states from delocalized 
wetting-layer states. In order to determine this approxi- 
mative energetic cutoff, which also limits the number of 
single-particle states considered as a basis for the FCI cal- 
culations later, we calculate e(kj_) for the polar (c-plane) 
and the nonpolar (m-plane) wetting-layer (without the 
QD) using Eq.([3|. We end up with 8 subbands in each 
case due to the sp 3 basis and the fact that the wetting- 
layer (WL) consists of two monolayers, as can be seen in 
Fig.Q and in Fig.([2| where the subbandstructures are 
depicted for their corresponding irreducible path. In the 
polar case, the initial sixfold symmetry of the Bravais- 
lattice is preserved while for the nonpolar wetting-layer 
the symmetry is reduced to a primitive rectangular one 
due to the fact that the structure has been rotated. The 
corresponding coordinates of the points of high symmetry 
are explicitely given in Tab.(|l]). To further illustrate this 
fact, the two-dimensional dispersion for the lowest con- 
duction subband is plotted in Fig.(fT|) for the polar and 
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TABLE I. Points of high symmetry for the polar (left) 
and nonpolar (right) wetting-layer with hexagonal (Cr >v ) and 
primitive rectangular (C^) symmetry respectively. Not given 
vector components ki are treated in real-space. 
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FIG. 1. Subbandstructure of the wetting-layer in the po- 
lar (c-plane) growth direction. The dispersion has hexagonal 
symmetry. As an inset the lowest two-dimensional conduction 
subband including the irreducible path is depicted. 



Fig. ([2]) for the nonpolar wetting-layer as an inset. The 
energetic cutoff for the electrons is now estimated to be 
the minimum of the conduction subband dispersions, i.e. 
min(e^^ =0 ) and the maximum of the valence subband 
ones, i.e. wax(e^ B =0 ) for the holes, respectively. This 
estimate is evaluated for both growth directions seper- 
ately and is justified, because if the QD is considered as 
only a small perturbation on an infinitely large wetting 
layer, these eigenenergies should in principle give the cor- 
rect energetic limit where delocalized wetting-layer wave 
functions should start to appear. In the polar (nonpolar) 
case we end up with an energetic cutoff of «1.821 (2.324) 
cV for the electrons and «0.429 (0.343) eV for the holes. 



B. Effective Confinement 

We start by presenting an excerpt of the calculated in- 
trinsic electrostatic potential energy in Fig.(3| according 
to Eq.( 12 ) obtained by the methods outlined in Sec. 3.B, 
as it is required to calculate the effective confinement 
potential. At first, there is a reduction of the maximum 




y s 

Irreducible wedge 



FIG. 2. Subbandstructure of the wetting-layer in the nonpo- 
lar (m-plane) growth direction. The dispersion has primitive 
rectangular symmetry due to a change in the shape of the 
unit cell. As an inset the lowest two-dimensional conduction 
subband including the irreducible path is depicted. 



0.15 




potential energy difference for the nonpolar orientation 



0.15 



FIG. 3. Intrinsic electrostatic potential energy in eV for the 
polar (c-plane) and nonpolar (m-plane) quantum dot orien- 
tation. Please note, that the actual numerical cell was much 
larger (f» 23.1 x 23.1 x 9.3 nm 3 ) as only the relevant excerpt 
is visualized. 



of about 100 meV from 488 meV to 384 meV compared 
to the polar case. This fact can be understood by the 
reduced spontaneous polarization part, since the surface 
area in the polar direction is also smaller. In each ori- 
entation, both fields tend to strongly seperate electron 
and hole wavefunctions into different directions. While 
in the polar case the electronic wavefunctions favor a 
localization in the top of the QD and the hole ones in 
the bottom near the wetting-layer, in the nonpolar ge- 
ometry, both carriers should prefer the opposite sides of 
the QD lens. The above argumentation assumes ground- 
state wave functions, but related to this assumption two 
questions arise: Does this also apply for excited states 
and, what is the role of many-body effects ? 

To address the first point, having calculated the in- 
trinsic electrostatic potential energy V p (R) = — e • (/>(R), 
it is now possible to determine an effective confinement 
potential by adding V^(R) as a diagonal contribution to 
Eq.(|3]). The effective confinement allows a first estimate 
of the energetic positions of the bound QD states in the 
heterostructure and provides more information on spatial 
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°0 10 20 30 40 50 60 

x-axis (a) 

FIG. 4. Side view of the effective confinement potential 
V e s(x) for the polar (c-plane) quantum dot for electrons (top) 
and holes (bottom). Each single line is a profile for different 
y-coordinates. The black solid line represents the electronic 
cutoff energy of ~ 1.821 eV while the dashed one corresponds 
to the hole cutoff of «0.429 eV. 




°0 5 10 15 20 25 30 35 40 
z-axis (c) 

FIG. 5. Side view of the effective confinement potential 
V e s(z) for the nonpolar (m-plane) quantum dot for electrons 
(top) and holes (bottom). Each single line is a profile for 
different y-coordinates. The black solid line represents the 
electronic cutoff energy of «2.324 eV while the dashed one 
corresponds to the hole cutoff of s^O.343 eV. 



localization perpendicular to the growth direction. 

Obtaining this field is achieved by calculating the elec- 
tronic dispersion according to Eq.([3| for every spatial 
in-plane position Rj^ of the full three-dimensional prob- 
lem, i.e. the quantum dot on the wetting-layer, but by 
using the actual material composition along Ry in real 
space and by adding the position dependant Vp(R|| ) as an 
additional diagonal contribution. This procedure corre- 
sponds to a solution of the full three-dimensional problem 
by replacing the potential variation perpendicular to the 
growth direction by a translational invariant potential, 
but solving the Hamiltonian for each position seperately. 
The effective confinement for electrons and holes in the 
QD can now be approximated by plotting the electronic 
and hole ground state energy of the WL problem as a 
function of position perpendicular to the growth direc- 
tion after exact diagonalization of Eq. |3| . 

Let us begin by discussing the polar (c-plane) quantum 
dot case for which the effective confinement potential is 
depicted in Fig.Q. One can expect the electronic wave- 
functions with low energy to be localized in the top of 
the quantum dot and the holes in the bottom respec- 
tively. For both electrons and holes the potential seems 
to be of flat and quite homogeneous nature and we ex- 
pect the ground states to be confined spatially in the QD. 
Excited states remain to be localized inside the QD and 
then gradually move to the wetting-layer. Remarkably, 
the effective electron potential is enhanced around the 
quantum dot for states with even higher energy than de- 
localized wetting-layer states, which is not the case for 
the holes as the electrostatic potential decays monotoni- 
cally. The spatial variation of V p (x,y), i.e. perpendicu- 
lar to the growth direction, is symmetric as it should be 



since the quantum dot has a sixfold rotational symmetry. 
Thus, there is no region which favors or supports local- 
ization of carriers due to the lack of xy-anisotropy in the 
effective confinement inside the quantum dot. As a re- 
sult, excited states can only gain overlap if they overcome 
the intrinsic fields along the z-direction. 

The results for the nonpolar geometry depicted in 
Fig. ([5]) show the modification of the carrier confinement 
due to the intrinsic electrostatic fields caused by the 
change in orientation. Generally, the bound hole states 
are localized in the lower part (z < zq) of the QD, while 
the electronic wave functions are pushed to the upper 
part (z > zq). Thus, electron and hole wavefunctions are 
separated for low energies and only excited states with a 
sufficient energy allow an increased overlap in the center 
of the QD. This is caused by the sloped carrier confine- 
ment which cleary shows the yz-anisotropy of V p (y, z), 
so that excited states can gain more overlap by mov- 
ing perpendicularly to the growth direction opposed to 
the polar case. Furthermore, one can observe that be- 
cause of this anisotropy the wetting-layer starts to be- 
come energetically more favorable than parts the QD 
if the energy is sufficiently high. The overall energetic 
position of the effective confinement potential V c ff(R) is 
shifted to a deeper position in the polar case with more 
strongly bound states. This behavior is expected as the 
field strength is enhanced in the polar orientation. 

As an additional benchmark for the energetic cutoffs 
calculated earlier, one can of course compare them to the 
present calculations including the intrinsic electrostatic 
potential. The energetic cutoffs in the polar case are 
plotted as the black solid line for the electrons («1.821 



eV) and as a dashed one for the holes («0.429 eV) in 
Fig.Q. In the case of the electrons, the approximative 
value obtained by neglecting the intrinsic fields and the 
quantum dot is by 6.7 meV lower than the value obtained 
with this more accurate calculation. In contrast, the hole 
cutoff is underestimated by about 39.1 meV which can be 
observed from Fig.Q as it does not approach the poten- 
tial in the limit of no QD. This fact indicates that the hole 
confinement is more sensitive to intrinsic fields compared 
to electrons, which is surprising, since one would expect a 
more remarkable deviation for both cases if the fields and 
quantum dot are neglected. The error introduced by this 
approximation only results in taking into account more 
single-particle hole states for the many-body calculation 
as actually to consider according to our approximative 
cutoff criteria. Note, that there is still some freedom 
in choosing the number of states, because the complete 
influence of the QD itself is not fully covered in these 
calculations using the factorization in Eq.(|3]). 

Furthermore, the energetic cutoff in the nonpolar ori- 
entation denoted by the black solid line for the electrons 
(«2.324 eV) and the dashed one for the holes («0.343 eV) 
is visualized in Fig. ([5]). In both cases the approximative 
value is by 1.2 meV lower than the value obtained by 
this more exact calculation and the shift is energetically 
symmetric for both carriers as expected. As can be seen, 
the approximative cutoff is a reasonable choice to distin- 
guish between quantum dot and wetting-layer states as 
it reproduces almost exactly the correct energetic limit 
of having no quantum dot and no intrinsic fields. 



C. Single-particle quantum dot states 

In the following we support the conclusions drawn in 
the previous analysis by taking an explicit look at some 
exemplarily chosen single-particle wave functions. Rep- 
resenting a general type of wave function in terms of spa- 
tial localization, we have plotted some manually selected 
ones for the polar and nonpolar orientation respectively 
in Fig. ^ and Fig. In the polar case, all states obey 
the sixfold rotational symmetry of the underlying lattice. 
Furthermore, the electronic states, from the ground-state 
up to a highly excited state, are localized inside the quan- 
tum dot and move gradually from the top to the bottom. 
In contrast, the degenerate hole ground-states are deeply 
localized in the wetting-layer, while excited states tend 
to be more extended around the QD but remain inside 
the WL. Only the highly excited hole states exhibit a 
significant probability density inside the QD. 

In contrast, the situation is completely different in the 
non-polar case as we can observe from Fig. At first, 
we are only left with a reflection as a symmetry opera- 
tion with respect to the xz-plane at y=0 since the intrin- 
sic electrostatic field with anisotropy in the z-direction 
spoils the former Ci v symmetry of the problem. Further- 
more, due to the sloped carrier confinement for electron 
and holes, their ground-state wave functions are strongly 
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FIG. 6. Visualization of selected single-particle states by 
plotting isosurfaces of equal probability density in the polar 
(c-plane) geometry. The quantum dot is indicated by the 
black line and the three lowest electrplotted from top and side 
view. In the upper part electronic wave functions are shown 
with increasing energy from left to right, i.e. the ground-state, 
an excited state in the middle and an highly excited state on 
the right. The same applies for the hole wave functions in the 
lower part. 



separated, even stronger than for the polar orientation 
due to the larger extent of the QD in the direction of 
the field. If we now increase the energy, excited states 
develop a finite probability density towards the center 
of the QD since the effective confinement now allows a 
localization with less separation. The higher excited elec- 
tron states have a finite probability distribution through- 
out the whole QD region, while the hole states seem to 
be restricted to the lower half of the QD. This can be 
understood by the fact that the effective hole mass is 
larger than the electron one, i.e. the holes are localized 
more strongly. Finally, we take a look at highly excited 
states which are supposed to feature localization inside 
the WL according to the discussion of the effective con- 
finement. Exactly this is the case for both electrons and 
holes, which still have a non vanishing probability density 
inside the QD, but their major density is located right 
above and below the QD in combination with a large 
distribution over the WL plane. 



D. Coulomb interaction vs. QCSE 

In this last section we extend the single-particle de- 
scription of the system to the many-body treatment out- 





FIG. 7. Visualization of some selected single-particle states 
by plotting isosurfaces of equal propability density in the non- 
polar (m-plane) geometry. The quantum dot is indicated by 
the black line and each state is plotted from side- and topview. 
In the left part electronic wavefunctions are shown with in- 
creasing energy from top to bottom, i.e. the ground-state, 
an excited state and finally a highly excited state. The same 
applies for the hole wave functions in the right part. 



lined earlier. The interplay between the QCSE and the 
Coulomb attraction manifests itself in the optical absorp- 
tion spectra. In the following, linear optical spectra in 
connection with the excitation of one electron-hole pair 
are determined as a function of the parameter 7 for the 
strenghth of the Coulomb interaction, Hx = Hq + jHc- 
The results of these calculations are depicted in Fig. ^ 
and Fig.(|9| and resemble the single-particle absorption 
spectra without Coulomb interaction in the case of 7 = 0, 
the true interaction strength for 7=1 and an unphysical 
artificial enhanced strength for 7 > 1. For visualization 
purposes we have plotted the logarithm of the absorption. 
For a system where Coulomb interaction only provides an 
energy renormalization, e.g. if it cannot compensate the 
intrinsic fields, and the electron and hole wave functions 
do not alter, we would expect that the energetic position 
of the emission line is reduced linearly with 7 if the inter- 
action is increased. At the point where Coulomb interac- 
tion becomes dominant over the intrinsic fields and allows 
electrons and holes to gain energy by reducing their sep- 
aration, i.e. they become spatially correlated, nonlinear 
effects should appear. Thus, the energetic position of the 
absorption lines should also be modified in a nonlinear 
way. We note, that a very small Lorentzian broadening 




0.5 1 1.5 2 2.5 
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FIG. 8. Absorption spectra of the polar quantum dot as a 
function of Coulomb interaction strength 7. 




interaction strength 



FIG. 9. Absorption spectra of the nonpolar quantum dot as 
a function of Coulomb interaction strength 7. 

has to be used for a proper analysis of the effect. We 
start by discussing the polar results in Fig.([8| and ob- 
serve that there are many emission lines with different 
slopes in the spectra due to different effective electron- 
hole pair distances. Nevertheless, the majority of these 
lines follow an almost linear trend for 7 6 [0,1]. At an 
artificially enhanced Coulomb interaction of 7 ~ 1.25 
several lines seem to split up into separate lines and non- 
linear behavior is clearly identified for 7 > 1.5. This 
is due to the fact that new excitonic states form in the 
system and the intrinsic fields are compensated by the 
Coulomb interaction. Please note that new bright lines 
start to appear in the spectra due to an increased interac- 
tion strength, which lifts degeneracies and mixes bright 
single-particle states into former dark excitonic states. 
In the nonpolar case depicted in Fig. ^ the situation is 
quite different, since most lines exhibit an identical lin- 
ear slope. As can be seen, in the large energy regime a 
splitting of lines appears at about 7 w 0.5 — 0.6 and the 
ground state transition clearly follows a nonlinear behav- 
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ior at 7 = 0.8. In addition there is a small indication of 
anticrossing of the absorption lines which also hints on 
strong energetic hybridization. Summarizing both cases, 
we note that the unmodified Coulomb interaction is able 
to compensate the intrinsic fields in the nonpolar geom- 
etry and 7 = 0.8 is the critical value according to our 
calculations. In contrast, only an unphysically enhanced 
interaction strength of 7 > 1.5 can compensate the QCSE 
in the polar geometry. To further illustrate this fact, we 
have visualized the corresponding many-body and single- 
particle ground state eigenfunctions according to Eq.( 16 ) 
for both geometries in Fig. ( 10 1 for the case of 7 = and 



7 = 1 including the corresponding absorption spectra rel- 
ative to the ground-state transition energy. In contrast, 
in this figure a lorentzian broadening of 10 meV for el- 
evated temperatures, Ref . (fTTj) . has been included which 
simulates the effect of all inelastic scattering processes 
not explicitly taken into account. 

As the Coulomb interaction is not strong enough to 
compensate the intrinsic electrostatic fields in the polar 
crystal orientation, there are only slight modifications of 
the excitonic ground-state wave function in comparison 
to the single-particle case. The main benefit is an energy 
gain due to the attractive electron-hole interaction. On 
the other hand in the nonpolar crystal orientation one 
can clearly observe that there is energy gain and also a 
strong spatial correlation of electrons and holes. As a re- 
sult, the many-body ground-state wave function is mod- 
ified drastically due to contributions of the formerly dis- 
cussed excited single-particle states with increased over- 
lap. Because the Coulomb interaction is strong enough 
for 7 > 0.8 to compensate the intrinsic fields, a new ener- 
getically more favorable ground-state is formed by mixing 
the excited singe-particle states. Thus, we end up with 
an enhanced absorption spectra in the nonpolar orienta- 
tion. Additionally, there is an overall energetic shift of 
the ground-state absorption peak in both crystal orienta- 
tions present. The corresponding emission spectra (not 
shown) have been published in Ref. (|3"rJ)) . 



lar wurtzite InN/GaN lens-shaped quantum dots. This 
was achieved by combining continuum elasticity the- 
ory and a tight-binding model to describe elastic and 
single-particle electronic properties of these nanostruc- 
tures properly. Before performing the full configuration 
interaction calculation for one clcctron-holc pair in or- 
der to solve the many-body problem, we calculated an 
effective confinement potential and discussed the single- 
particle states concerning their spatial localization inside 
the heterostructure on that basis. As a partial result, the 
sloped carrier confinement in the nonpolar geometry led 
to a different localization behaviour of excited states. In 
detail, electrons and holes tend to be less separated for 
excited states and surprisingly even higher excited states 
move spatially to the wetting-layer. The treatment of the 
interacting many-body problem of one exciton revealed, 
that the Coulomb interaction can compensate the intrin- 
sic electrostatic fields and the critical relative strength 
is 7 = 0.8 in the nonpolar quantum dot. On the other 
hand, the Coulomb interaction is not strong enough to 
compete with the polarization fields in the polar quan- 
tum dot geometry, since only an artificially enhanced 
relative strength of 7 = 1.5 compensation took place. 
These insights were achieved identifying a nonlinear ef- 
fect in the energies of optical transitions as a function of 
Coulomb interaction strength. Furthermore, the nonpo- 
lar ground-state absorption peak exhibits a quantitative 
enhancement due to Coulomb interaction and not only 
an energetic shift as in the case in the polar orienta- 
tion. According to these findings, practical consequences 
for building blocks are related to higher carrier-densities, 
because in that regime one could expect multiexcitonic 
complexes in nonpolar structures to compensate the in- 
trinsic fields much faster than in their polar counterparts. 
This could lead to an enhancement of optical recombina- 
tion processes in that regime. 



V. SUMMARY 



In this work we investigated the interplay between 
Coulomb interaction and the QCSE in polar and nonpo- 
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